Bloch-type Domain Walls in Rhomb ohedral BaTiOs 
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Ferroelectric domain walls (FDWs) are usually considered to be of Ising type, but there have 
been suggestions in recent years that Bloch-type FDWs are also possible in some cases, e.g., in 
the rhombohedral phase of BaTiOs. The mechanically compatible and electrically neutral FDWs in 
rhombohedral BaTiOs are of 71°, 109°, and 180° type. We have investigated these FDWs based both 
on first-principles calculations and on a Ginzburg-Landau-Devonshire (GLD) model [P. Marton, I. 
Rychetsky, and J. Hlinka, Phys. Rev. B 81, 144125 (2010)]. The results from both approaches 
confirm the Ising nature of the 71° FDW and the Bloch nature of the 180° FDW, and predict both 
Ising-type and Bloch-type FDWs are possible for the 109° case. Considering the relatively small 
rhombohedral strain in BaTiOs, the competition between the energies of Bloch and Ising FDWs can 
be discussed in terms of a picture in which a Bloch wall is regarded as being composed of a pair of 
smaller-angle Ising ones. A reduction by 40% in the parameters describing the gradient term in the 
GLD model brings it into better agreement with the first-principles results for detailed properties 
such as the energies and widths of the FDWs. 

PACS numbers: 77.80.Dj,77.22.Ej,77.84.-s 



I. INTRODUCTION 



Ferroelectrics find many industrial and commercial ap- 
plications, such as in high-dielectric constant capacitors, 
ferroelectric thin-film memories, piezoelectric transduc- 
ers, nonlinear optical devices, and switches. The perfor- 
mance of many kinds of ferroelectric devices is affected 
by the ferroelectric domain structure and the properties 
of the domain boundaries. For example, a recent experi- 
ment has shown that the observed dielectric permittivity 
of a BaTiOs single crystal in the rhombohedral phase 
varies depending on the domain structures induced by 
pre- treatments at higher temperatures. 2 Such influences 
on the mechanical and electrical properties of devices has 
motivated theoretical and experimental work directed to- 
ward obtaining a better understanding of ferroelectric 
domain structures. 

Ferroelectric domain walls (FDWs) are usually con- 
sidered to be of Ising type, in which Py, the projection 
of the polarization vector onto the plane of the domain 
wall, simply reverses itself by passing through zero along 
a high-symmetry path as one scans through the domain 
wall. Ising FDWs tend to be favored because ferro- 
electrics are generally strongly electrostrictive, so that a 
rotation of Pn away from this high-symmetry path would 
entail a significant elastic energy cost. In contrast, the 
spontaneous magnetostriction which couples the magne- 
tization and lattice strain in ferromagnetic materials is 
typically much weaker. As a result, magnetic domain 
walls are usually much wider, on the order of microns, 
and the magnetization vector can rotate away from the 
high-symmetry path. The domain wall is denoted as a 
Bloch or Neel wall depending on whether this rotation 
occurs in a plane parallel or normal to the domain wall, 



respectively. 

In recent years, however, there have been some the- 
oretical predictions of the presence of Bloch and even 
Neel components in some ferroelectric materials and 
heterostructures. 3 5 In particular, it has been predicted, 
in the framework of a phenomenological Ginzburg- 
Landau-Devonshire (GLD) model, that the 180° FDWs 
in rhombohedral BaTiOs should be of Bloch type. 6 This 
work has motivated us to test whether this behavior is 
also reproduced by first-principles density-functional cal- 
culations on BaTiOs. 

Since the discovery of ferroelectricity in this mate- 
rial in 1954, 1 BaTiOs has been very widely studied and 
has emerged as a kind of prototypical ferroelectric com- 
pound. It undergoes a sequence of phase transitions from 
a high-temperature paraelectric cubic phase to ferroelec- 
tric tetragonal, orthorhombic, and finally rhombohedral 
phases as the temperature is reduced. Here we are in- 
terested in the zero-temperature rhombohedral phase, in 
which the spontaneous polarization prefers to lie in eight 
energetically equivalent directions, as shown by the ar- 
rows in Fig. 1. The figure also shows the possible rota- 
tion angles between the spontaneous polarization direc- 
tions on the two sides of the domain wall (relative to the 
arrow marked as '0'). In the low- temperature rhombohe- 
dral phase of BaTiOa, the FDWs are therefore of three 
types: R71°, R109°, and R180°. (The 'R' denotes a 
FDW in the rhombohedral phase, following the notation 
of Ref. 6.) Taking into account the constraints of electri- 
cal neutrality and mechanical compatibility, the plane of 
the FDW is normal to the sum of the two polarization 
vectors for the R71° and R109° cases, while for the 180° 
case it can be either {211} or {110}. 

In this work the R71°, R109°, and R180°{110} FDWs 



FIG. 1. (Color online) The directions of the symmetry- 
allowed spontaneous polarizations in rhombohedral BaTiOs. 
Angles are relative to the reference direction labeled as 0°. 



in BaTiOs are investigated using first-principles calcula- 
tions in the context of density-functional theory. While 
first-principles calculations have been successfully ap- 
plied to study FDWs in many other cases, 3,5,7,8 we are 
not aware of any previous application to the case of the 
rhombohedral phase of BaTiOs. Our calculations con- 
firm the Bloch nature of the R180°{110} FDW and the 
Ising nature of the R71° FDW, in agreement with the 
predictions of the GLD model. The energy difference be- 
tween the Bloch R109° and Ising R109° FDW is very 
small, which suggests that both types of the domain wall 
are possible. The Bloch walls can be regarded as being 
composed of two Ising walls of smaller angles. A compar- 
ison of the sum of the energies of these constituting Ising 
walls with that of an Ising-type solution can explain why 
the polarization vector picks a Bloch-type path in some 
FDWs. A quantitative comparison of our first-principles 
results with those of the GLD model suggests that a 40% 
reduction in the size of a gradient term in the GLD model 
is needed to bring the two theories into good agreement. 

The manuscript is organized as follows. In Sec. Ill we 
describe the geometry of each of the FDWs to be stud- 
ied. We also review the first-principles and GLD model 
approaches which are used to study the FDWs, and give 
the details of the methods used for the first-principles 
calculations. The results from the first-principles calcu- 
lations and their comparison to the GLD model are then 
described in Sec. IV. In Sec. V we discuss the competi- 
tion between Ising and Bloch configurations in terms of 
energy considerations, and in Sec. VI we briefly summa- 
rize and discuss future prospects. 



II. DOMAIN WALL AND SUPERCELL 
GEOMETRIES 

The mechanically compatible and electrically neutral 
FDWs investigated in this paper are shown in the left 
column of Fig. 2, where the arrows indicate the orienta- 
tion of the polarization vectors P(— oo) and P(oo) on the 
two sides of the domain wall. In the right column, the 
symmetry- adapted coordinate system (r, s,t) is shown for 
each of these walls. The unit vector normal to the wall 
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FIG. 2. (Color online) The R71°, R109° and R180°{110} 
FDWs in BaTiOs. The arrows in the figures at left show 
the directions of the polarization vectors on the two sides 
of the FDW, while those on the right depict the associated 
symmetry- adapted coordinate system (r, s,t). 



is denoted by s. The second unit vector r is chosen to be 
parallel to P(oo) — P(— oo), the difference between the 
spontaneous polarizations on the two sides of the wall; 
electrically neutrality implies that this is normal to s. 
The third basis vector is defined as t = r x s. 

For the application of the GLD continuum approach, 
atomistic details are not important, and specific atom- 
istic geometries do not have to be considered. This is 
obviously not the case, however, for the first-principles 
calculations. These are set up by considering a super cell 
that is extended along the direction s normal to the wall, 
keeping minimal dimensions in the orthogonal directions. 
Ideally we would prefer a supercell containing only a sin- 
gle domain wall, but this is incompatible with periodic 
boundary conditions. Thus, we use supercells contain- 
ing two equivalent FDWs, where we enforce this equiv- 
alence by imposing a two-fold screw symmetry (two-fold 
rotation about s followed by a half superlattice-vector 
translation along the stacking direction s). If these walls 
are well separated, their effect on each other should be 
small and the physical properties of the walls should be 
unaffected. 
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FIG. 3. (Color online) 10-atom rotated building block that is 
stacked along s to construct supercells for studying R71° and 
R180°{li0} FDWs, which lie normal to s. 
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FIG. 4. (Color online) A supercell with four rotated 10-atom 
units stacked in the s direction. The centers of the FDWs are 
shown by the orange planes. 



We first construct a reference paraelectric supercell by 
identifying a minimal building block having lattice vec- 
tors parallel to r, s and t, and repeating this block N 
times along the stacking direction s. An example of 
such a building block, used for the R71° and R180°{110} 
cases, is shown in Fig. 3, and an example of a supercell 
built from it is shown in Fig. 4. An initial configuration 
for an Ising FDW is then chosen by shifting the coor- 
dinates of the oxygen atoms along the P(— oo) direction 
in the first half of the supercell, and along P(oo) in the 
second half. For the R71° and R109° cases this results in 
a configuration with a mirror symmetry relating the — t 
and t directions, and since this symmetry is preserved by 
the subsequent relaxation of atomic coordinates, the re- 
sulting relaxed configuration is guaranteed to be of Ising 
type. To initialize a calculation on a Bloch FDW, we 
also add oxygen displacement components along t in one 
FDW and along — t in the other (still preserving the screw 
symmetry) , so that effectively the displacement vector in 
the r-t plane is rotated from the r to the t direction 
across the wall before pointing to -r on the other side 
of the wall. In both the Ising and Bloch cases, the s 
component is left unchanged in the initial configuration, 
although of course it may relax later. 

We then relax the atomic coordinates until all of the 
forces fall below a chosen threshold. While doing this, 
we constrain the two in-plane lattice vectors (i.e., in 
the plane of the FDW) to remain fixed, consistent with 
the relaxed strain state of a single-domain rhombohe- 
dral crystal. We do this because the physical system we 
are trying to model is an isolated FDW between very 
thick domains, in which case the bulk elasticity domi- 
nates over the interface and fixes the in-plane strain. We 
let the third (long) superlattice vector relax along with 
the atomic coordinates during the minimization. Finally, 
the polarization profile is calculated from the pattern 
of displacements and the calculated dynamical effective 
charges as described in Sec. Ill A below. 

For all the investigated FDWs, we can always find 
an Ising-type FDW solution if we impose an appropri- 
ate symmetry constraint. To test whether this solu- 
tion is locally stable, we add small symmetry- lowering 
atomic displacements, and check whether it relaxes back 
to the Ising solution. We next try starting from a Bloch- 
type configuration, with substantial distortions similar 
to those that would be present in well-defined domains 



of rhombohedral phase. Again, if this structure relaxes 
back to the Ising one, then we conclude that no Bloch 
FDW was found, and the Ising solution is stable. If we 
find instead that the calculation converges to a Bloch-like 
solution, we compare its energy with that of the Ising so- 
lution (if locally stable) to determine which is the global 
solution. 

The R71° wall has its normal s in the [110] direction, 
so that it lies parallel to the diagonal plane in the prim- 
itive cell as shown in Fig. 2. Therefore, we consider the 
building block of Fig. 3, which is a 10-atom x x 1 
cell obtained by rotating by 45° around the z axis with 
respect to the parent cubic unit cell. In this rotated cell, 
the FDW lies in the (100) plane. The simulation supercell 
is constructed by stacking these units in the s direction, 
as shown in Fig. 4. The initial coordinates for the Ising 
R71° FDW simulation are then obtained by shifting the 
oxygen atoms by about 0.1 A along [101] in half the su- 
percell, and along [101] in the other half, in the rotated 
coordinate system. 

The R109° FDW lies in the (100) plane, s = x, and 
a supercell can easily be made by stacking the primi- 
tive 5-atom rhombohedral cells in this direction. The 
initial Ising configuration is then set by displacing the 
oxygen atoms from the equilibrium positions along the 
[111] direction on one side of the FDW, and along the 
[111] direction on the other side. 

The R180°{110} wall is again parallel to the diagonal 
plane in the primitive cell, as for the R71° domain wall. 
So in this case as well, a supercell is made by stacking the 
rotated 10-atom units shown in Fig. 3 in the s direction. 
The initial configuration for the R180°{110} Ising FDW 
is obtained by shifting the oxygen atoms from their equi- 
librium positions along the [011] direction on one side of 
the FDW and along the opposite direction on the other 
side, in the rotated coordinate system. 

For the Ising R180°{110} case there is also the possi- 
bility of imposing a higher symmetry by insisting that the 
FDW lie exactly in a Ba-Ti-0 plane, or exactly in an 0-0 
plane, which can be accomplished by adopting an inver- 
sion center through one of the atoms lying in the FDW. 
This is not possible for the R71° and R109° FDWs, since 
the presence of a global P s component makes the "front" 
and "back" sides of the FDW distinguishable and rules 
out the presence of an inversion symmetry. 



III. COMPUTATIONAL APPROACH 
A. First-principles calculations 

The calculations are done using the ABINIT imple- 
mentation of density functional theory 9 within the local- 
density approximation (LDA) using the Perdew-Zunger 
exchange-correlation functional. 10 Ultrasoft pseudopo- 
tentials, 11 in which semicore s and p states are included 
in the valence for Ba and Ti, were converted for use as 
projector augmented- wave (PAW) potentials 12 using the 
USPP2ABINIT package. The plane wave cut-off and the 
energy cut-off for the fine FFT grid are set to 25 Ha 
and 40 Ha respectively. The tolerance on the difference 
of forces in successive iterations in a self-consistent-field 
(SCF) cycle is set to 5.0 x 10~ 10 hartree/Bohr, which 
reached twice successively, causes one SCF cycle to stop 
and ions to be moved. The structural optimizations 
are done using the Broyden-Fletcher-Goldfarb-Shanno 
minimization. 13 

The supercells employed for studying R71° and R180° 
FDWs are of dimensions y/2Na x y/2a x a, where a is 
the primitive lattice constant and N is the number of 10- 
atom units stacked in the s-direction. Therefore, a 1 x 
4x6 Monkhorst-Pack 15 fc-mesh is chosen for simulating 
these domain walls. The supercell for the R109° FDW 
is made by stacking the 5-atom primitive cells in one 
direction, so a 1 x 6 x 6 Monkhorst-Pack 15 /c-mesh seems 
a proper choice for simulating this FDW. 

We also compute the dynamical effective charges 16 Z* 
in bulk paraelectric cubic BaTiOs. The dynamical charge 
tensor Z* a a of a given atom measures the dipole induced 
along f3 by a displacement of atom i along a. In many 
oxides including BaTiOa, these charges are quite differ- 
ent from the formal ionic charges. The Z* tensors are 
computed by finite differences, i.e., by making small dis- 
placements and calculating the resulting change in Berry- 
phase polarization. 17 

These dynamical charges are then used as an ingredient 
in an algorithm 3 by which we map out the polarization 
profiles in FDW-containing supercells, as follows. The 
polarization is only changing along the stacking direction 
and is constant in the planes normal to this direction. So 
first the contribution from each layer to the polarization 
in direction a arising from displacements of atoms j in 
direction j3 is calculated as 

P,jei 

where I is a layer index. In the supercell used for study- 
ing the R109° FDW these are Ba-0 and Ti-O-0 layers, 
while for the R71° and R180°{110} FDWs built from the 
rotated 10-atom units I refers to Ba-Ti-0 and 0-0 lay- 
ers. If we break the supercell into smaller cells centered 
on these layers, we can assign a local polarization to cell 
I by counting its own contribution and half that of each 



neighbor, i.e., 

Pi° = ^(^- 1) +^ ) + ^ +1) ), (2) 
where Q is the volume of the cell. 



B. The GLD model 

We review the Ginzburg-Landau-Devonshire model 
used in Ref. 18, which is again used here to model the 
FDW properties and compare with the first-principles re- 
sults. The excess free energy F relative to the reference 
cubic paraelectric state is expressed in terms of polariza- 
tion and strain fields as 

F[{P i ,P iJ ,e iJ }] = J f(r)dr, (3) 

where / is the GLD free-energy density which is taken 
to be a function of the polarization components P^, their 
spatial derivatives Pij = dPi/dxj, and strain compo- 
nents Cij. In particular, / is expressed in terms of Lan- 
dau, elastic, electrostriction, and gradient terms: 

/ = f { L } {Pi} + fdeij} + f q {Pi, e« } + f G {Pi,j}- (4) 

In Ref. 18, explicit forms were given for each of the terms 
in this expression, parameter values were estimated from 
the bulk single-crystal properties of BaTiOa, and the 
GLD model was used to investigate domain-wall prop- 
erties. 

Here we are especially concerned with the gradient or 
Ginzburg terms in the free-energy expansion, which take 
the form 

/ G =^G 11 (P 1 2 1 + P 2 2 2 + J P 3 2 3 ) 

+ Gi, 2 (Pi,iP 2> 2 + ^2,2^3,3 + Pl,lPz,z) 

1 ( 5 ) 

+ -G44 [(Pl,2 + P2,l? + (P 2 ,3 + P 3 , 2 ) 2 
+ (P3,1+Pl,3) 2 ] • 

As discussed in Ref. 18, considerable caution was required 
in extracting the G coefficients from inelastic neutron 
scattering experiments, and the remaining uncertainties 
are significant. The G tensor has an important effect on 
the widths and energies of the FDWs, so that the uncer- 
tainties in the values of these coefficients is a limiting fac- 
tor in determining the properties of the FDWs. The orig- 
inal parameters of Ref. 6 describing the gradient terms 
are Gn = 51xl0- n Jm 3 C- 2 , G 12 = -2x lO" 11 Jm 3 C- 2 , 
and G 44 = 2 x 10- 11 Jm 3 C- 2 . 

In order to establish a better estimate of these coeffi- 
cients, the GLD model is employed below to recalculate 
the polarization profiles of R71°, R109°, and R180°{110} 
FDWs at zero temperature using modified G coefficients, 
and the results are compared with first-principles ones. 
In order to facilitate the comparisons, these GLD model 
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TABLE I. The experimental and theoretical values of the dy- 
namical effective charges of the Ba, Ti and atoms in BaTiOs 
(in units of the charge quantum e). 
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Z* (LDA ) 


2.75 


7.18 


-1.86 


-5.65 



calculations have been performed on the identical geome- 
tries as in the first-principles calculations. That is, we 
impose periodic boundary conditions corresponding to 
the size of the first-principles supercells, impose the same 
two- fold screw symmetry as was used there, and specify 
the strain state to be consistent with infinite domains. 



Q_ 

CL 



(a) 

5 

0- 


X 


5 




n 


X 



Q_ 



CL 
CL 
CL 



-1.0 -0.5 0.0 0.5 1.0 




Domain-wall width 



-1.0 -0.5 0.0 0.5 1.0 
Pr/P^ 



(b) 








/ 


— Pr 

~ Pt 
~ Ps 


V / 




\ / 


D12345678C 
s [nm] 


(d) 










— Pr 
~ Pt 
* * Ps 








2 4 6 8 

s [nm] 



In the case of a domain wall whose polarization profile 
can be fit to a hyperbolic tangent, p(x) = po tanh(#/£), 
a common definition of the width is w = 2£ . However, 
some of the domain walls to be studied here have un- 
usual polarization profiles that do not resemble a single 
hyperbolic tangent at all. To accommodate such cases, 
we define w as the width of the region within which 
\p(x)\/po < tanh(l) = 0.762. This definition has the 
advantages of being globally reasonably and of reducing 
to the conventional definition above when the FDW does 
resemble a hyperbolic tangent. We adopt this definition 
throughout the remainder of this work. 



IV. RESULTS 

BaTiOs has been studied extensively both experimen- 
tally and theoretically. Our first-principles computed val- 
ues of 3.95 A for the lattice constant and 89.93° for the 
rhombohedral angle can be compared with experimental 
values 19 of 4.00 A and 89.87°, respectively. Our results 
are consistent with the experience that the LDA typically 
gives slightly underestimated values for the unit-cell vol- 
ume and ferroelectric distortion. 20 

We have also computed the values of the dynamical 
charge tensors for cubic paraelectric BaTiOs, to be used 
as an ingredient in the algorithm for computing polar- 
ization profiles as described in Sec. Ill A. The calculated 
values of Z* for Ba, Ti and O in BaTiOs are compared 
to the experimental values as well as ionic charges in Ta- 
ble I. The values of Zo || and Zo JL refer to the Z* of the 
oxygen ion when it is displaced along the Ti-0 direction 
or perpendicular to it, respectively. 

A Berry-phase calculation of the polarization 17 using 
the relaxed atomic positions in rhombohedral BaTiOs 
yields a value of 30/iC/cm 2 . The polarization calcula- 
tions using the atomic displacements and the theoret- 
ical Z* values yields 31/iC/cm 2 , which is slightly un- 



FIG. 5. (Color online) Polarization profiles of Ising R71° 
FDWs in rhombohedral BaTiOs as calculated using (a-b) 
first-principles calculations, and (c-d) the GLD model. Pan- 
els (a) and (c): parametric plots showing polarization values 
in the P r —Pt plane. Crosses indicate local energy minima as- 
sociated with homogeneous rhombohedral domains. Panels 
(b) and (d): polarization components as a function position 
s along the supercell direction (two supercells containing four 
FDWs are shown for clarity). 



derestimated compared to the experimental value 19 of 
33.5/iC/cm 2 as expected. The GLD model, on the other 
hand, yields a value of 38/iC/cm 2 . 

In the following subsections, we report the results 
of our supercell calculations for the R71°, R109°, and 
R180°{110} ferroelectric domain walls. 



A. The R71° domain wall 

We have investigated the R71° FDW in rhombohedral 
BaTiOs by carrying out first-principle calculations on 
an 80-atom supercell made by stacking 10-atom rotated 
units as described in Sec. II. In this case we only found 
an Ising FDW; perturbing this by adding symmetry- 
lowering components only led back to the Ising structure 
upon further relaxation. The polarization profile com- 
puted from the relaxed Ising structure using Eq. (2) is 
displayed in two ways in Fig. 5 (a-b). The left panel, 
Fig. 5(a), shows the Pr^ and P^p values for each layer 

/, while Fig. 5(b) shows plots of Pr \ pP and P^p as a 
function of position s while scanning through a sequence 
of four domain walls (two entire supercells). As is clear 
from this figure, the P t component remains zero every- 
where in the supercell, which clearly indicates the Ising 
nature of this FDW. The P r value reverses quite sud- 
denly and attains a value very close to its saturation bulk 
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FIG. 6. (Color online) Polarization profiles for Ising R109° 
FDWs in BaTiOs using (a-b) first-principles calculations and 
(c-d) the GLD model. Details are as in Fig. 5. 
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FIG. 7. (Color online) Polarization profiles for Bloch R109° 
FDWs in BaTiOs using (a-b) first-principles calculations and 
(c-d) the GLD model. Details are as in Fig. 5. 



value deep inside each domain, indicating a rather narrow 
FDW width. The P s value is almost exactly constant; as 
we shall see, this is true of all FDWs in this study, con- 
sistent with the expectation that inhomogeneities in P s 
would result in bound charge which in turn would involve 
an extra Coulomb energy cost. 

The above results are in good qualitative agreement 
with those of the GLD model, but the GLD domain- 
wall width of 0.58 nm is almost twice that of the first- 
principles prediction of 0.33 nm. This suggests the need 
for a reduction of the zero-temperature values of the gra- 
dient terms in the GLD model to bring it into better 
agreement with the first-principles results. This reduc- 
tion of the gradient terms will be even more important 
in the R180°{110} case, to be discussed in Sec. IV C. As 
explained there, we have chosen to reduce all three of the 
gradient coefficients in Eq. (5) by 40% in order to arrive 
at an improved GLD model. 

The results computed for the R71° FDW using this 
modified GLD model are presented in Fig. 5 (c-d) using 
the same plotting conventions as for the first-principles 
results in Fig. 5 (a-b). The FDW width is now 0.37 nm. 
Furthermore, the GLD domain- wall energy is reduced 
from 5.0 to 3.2mJ/m 2 , to be compared with the first- 
principles value of 3.8mJ/m 2 . While these numerical 
values should be interpreted reservedly in view of the un- 
certainties in both theories, it is clear that the GLD the- 
ory is in better agreement with the first-principles theory 
after the reduction of the strength of the gradient term. 



B. The R109° domain wall 

According to the GLD model calculations of Marton 
et a/., 6 both Ising and Bloch solutions are possible for 
the R109° FDW in BaTi0 3 . Our first-principles results 
confirm this picture. 

Starting first with the Ising case, Fig. 6 (a-b) shows the 
polarization profile for a 50-atom supercell in which the 
atomic positions have been relaxed from an initial config- 
uration with two Ising-type R109° FDWs. The domain 
wall is again fairly narrow, though not quite as narrow 
as in the R71° case. The energy and width of this Ising 
R109° FDW are calculated from first principles to be 
0.36 nm and 11.1 mJ/m 2 , respectively. The correspond- 
ing GLD results using the reduced gradient term, shown 
in Fig. 6 (c-d), are clearly in good qualitative agreement. 

When the atomic positions are relaxed from an ap- 
propriately distorted initial configuration, a Bloch-type 
solution for this wall is found. The first-principles polar- 
ization profiles computed for the Bloch-type R109° FDW 
in a 50-atom supercell are shown in Fig. 7(a-b). To a 
first approximation, this Bloch FDW can be regarded 
as a composition of two 71° FDWs in close proximity. 
In part for this reason, the Bloch-type FDW is clearly 
broader than the Ising one. However, the energies of the 
Bloch and Ising solutions are found to be almost identi- 
cal, with the Ising one being only ^2-3% lower in energy. 
If we extrapolate to larger separations between FDWs 
we might expect the Bloch energy to fall more than the 
Ising one, because of the larger FDW width in the Bloch 
case. This suggests that both types of R109° FDWs 
have similar energies and that both might be found in 
rhombohedral BaTiC>3 crystals. The reason for existence 
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FIG. 8. (Color online) Polarization profiles for Ising 
R180°{110} FDWs in BaTiOs using first-principles calcula- 
tions. Details are as in Fig. 5. 



FIG. 9. (Color online) Polarization profiles for Bloch 
R180°{110} FDWs in BaTiOs using (a-b) first-principles cal- 
culations and (c-d) the GLD model. Details are as in Fig. 5. 



of both Ising and Bloch solutions for the R109° FDW 
is discussed in Sec. V. The modified GLD model again 
gives good qualitative and semiquantitative agreement 
with the first-principles results for the case of the Bloch 
R109° FDW, as shown in Fig. 7(c-d). 



C. The R180°{110} domain wall 

A major result of the GLD study of Marton et a/., 6 
was the prediction that the lowest-energy FDW for the 
R180°{110} case in BaTi0 3 is of Bloch type, with an 
energy lying 10% lower than that of its Ising counter- 
part. Verifying this result from first principles is com- 
putationally more challenging than studying the R71° or 
R109° FDWs. The R180°{110} FDW has the biggest 
rotational angle and is hence the widest of the three in- 
vestigated FDWs. Moreover, the Bloch- type FDWs are 
generally much broader than the Ising- type ones. Prelim- 
inary first-principles calculations showed that supercells 
smaller than 80 atoms are too small to accommodate 
two R180° FDWs; initial 60-atom Bloch-wall supercells 
relaxed to unreasonable configurations. We have there- 
fore carried out our calculations on an 80-atom supercell. 
While the polarization does not quite have room to reach 
its saturation value between neighboring FDWs, at least 
we obtain a stable solution that can reveal the Ising or 
Bloch nature of this FDW. 

As mentioned in the last paragraph of Sec. II, it is pos- 
sible to enforce an Ising-like geometry in the R180°{110} 
case by imposing an initial inversion symmetry about 
an atom in the center of the FDW and preserving 
this symmetry during relaxation. When we apply our 



first-principles calculations including this symmetry con- 
straint, we arrive at a configuration like that shown in 
Fig. 8 (a-b), which is Ising-like in the sense that P r = 
P t = in the center of the wall. However, P t has sub- 
stantial excursions away from zero, with the polarization 
path following an S-like curve in (P r ,P £ ) space, as can 
be seen clearly in Fig. 8(a). 

However, we find that if we do not impose this spe- 
cial symmetry, most initial conditions relax to the Bloch 
configuration shown in Fig. 9 (a-b). Even if we start from 
an Ising-like configuration and break the symmetry only 
slightly, we find that the simulation will eventually relax 
to the Bloch configuration. It is clear from Fig. 9 (a-b) 
that both P t and P r components are strongly non-zero, 
and the rotation of the polarization as one progresses 
through the domain wall, which is the characteristic fea- 
ture of a Bloch-type wall, is clearly visible. It is also 
obvious that this Bloch FDW resembles an adjacent pair 
of 71° and 109° FDWs. We also find that this Bloch wall 
has a significantly lower energy than that of the Ising- 
like wall of Fig. 8, supporting the conclusion that the 
Bloch solution is the global minimum for the case of the 
R180°{110} FDW. 

We have also carried out corresponding simulations of 
the Ising-like and Bloch configurations of the R180°{110} 
FDW using the GLD model. If this is done using the full 
strength of the gradient term under periodic boundary 
conditions, we find a stable Bloch solution only when 
the centers of the domain walls can be at least 3.4 nm 
apart. This is equivalent to using a supercell of 120 
atoms in the first-principles calculations, and imposing 
the rhombohedral epitaxial strain does not change this 
result very much. On the other hand, we find that if the 
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TABLE II. (Color online) Summary of computed energies and 
widths of R71°, R109°, and R180°{110} FDWs in BaTi0 3 . 
LDA indicates first-principles results; GLD and GLD-r refer 
to the Ginzburg-Landau-Devonshire model with original and 
reduced gradient term, respectively. 



Wall width (nm) Wall energy (mJ/m 2 ) 
LDA GLD GLD-r LDA GLD GLD-r 



Ising R71° 


0.33 


0.58 


0.37 


3.8 


5.0 


3.2 


Ising R109° 


0.36 


0.54 


0.34 


11.1 


10.6 


6.8 


Bloch R109° 


1.01 


1.10 


0.72 


11.2 


10.2 


6.2 


Ising R180° 


0.51 




0.83 


26.6 




27.2 


Bloch R180° 


1.38 




1.40 


24.0 




27.6 



FIG. 10. (color online) The path swept by the tip of the po- 
larization vector while rotating from Pi to P3 when crossing 
through hypothetical Bloch-type FDWs for the R71°, R109°, 
and R180°{110} geometries. For the R71° case the path 
(which rotates in a plane parallel to the FDW) is shown in 
red to indicate that it is unfavorable. 



R71° R109° R180°[110] 




gradient term is reduced by about 40%, then stable solu- 
tions become possible under the same periodic boundary 
condition as in the 80-atom supercell. The polarization 
profiles calculated for the R180°{110} FDWs using the 
GLD model with reduced gradient term, with the same 
domain-wall distance and elastic boundary conditions as 
in the first-priciples calculations, are shown in Figs. 8(c- 
d) and 9(c-d) for the Ising- like and Bloch cases respec- 
tively. The GLD results are clearly now in good qualita- 
tive agreement with the first-principles calculations. 

The computed energies and widths of all of the FDWs 
are collected and presented in Table II. It is evident 
that the first-principles and GLD results are in broad 
agreement. As expected, the GLD model with reduced 
gradient term yields narrower walls and lower wall ener- 
gies, yielding improved agreement for the Ising R71° and 
R109° cases, but somewhat overshooting for the Bloch 
R109° case. 

For the R180° FDWs, the values given in Table II 
ought not be taken too seriously because the repeat dis- 
tance of the supercell is rather short compared to the 
FDW width. In fact, we did not succeed in finding sta- 
ble FDW solutions with the original GLD model. The 
present 80-atom supercell is large enough to give sta- 
ble solutions in both the LDA and reduced-GLD cal- 
culations, but their properties are undoubtedly not yet 
converged with supercell size. The FDW energy can be 
expected to fall with supercell size, so the energies in 
Table II should be taken as upper bounds. The "crowd- 
ing" of the Bloch R180° walls appears to be more serious 
than for the Ising ones, which can explain why the re- 
duced GLD model predicts a slightly higher energy for 
the Bloch compared to the Ising FDW. When the GLD 
calculations are repeated for larger FDW separations, 
they clearly predict that the Bloch configuration is lower 
in energy. 6 The first-principles calculations already pre- 
dict the Bloch wall to be lower in energy for the 80-atom 
supercell, and this trend would only be strengthened if 
we could afford to repeat the calculations at increasingly 
larger FDW separations such that well-developed rhom- 
bohedral domains could form between domains. 



V. DISCUSSION 

Among the three investigated FDWs, the polarization 
vector rotates by the smallest angle in the R71° FDW, 
and by the biggest angle in the R180°{110} FDW, so it 
is not surprising that the former has the smallest energy 
and the latter the biggest, as summarized in Table II. We 
can also propose a simple explanation for the fact that 
the first-principles calculations and the GLD model re- 
sults predict an Ising nature for the R71° FDW, a Bloch 
nature for the R180°{110} FDW, and a very small en- 
ergy difference between Ising and Bloch solutions for the 
R109° FDW. The path the polarization vector would take 
in rotating from one side of the wall to the other in a 
Bloch-type solution is shown for each of these FDWs in 
Fig. 10. In a hypothetical Bloch R71° FDW the polar- 
ization vector would pass close to the center of one of 
the adjacent faces, which corresponds to a tetragonal po- 
larization state and is not energetically favorable in the 
rhombohedral phase. 

On the other hand, the Bloch 109° FDW can be con- 
sidered as a combination of two R71° FDWs. As can be 
seen from Table II, the total energy of two Ising R71° 
FDWs is rather close to the energy of one Ising 109° 
FDW, from both the LDA and GLD calculations. By 
way of a caveat, we point out that such a comparison 
may be overly simplistic because the Ising walls com- 
prising the Bloch 109° FDW experience a foreign strain 
environment and do not conform to the plane of mechan- 
ical compatibility of a true R71° FDW. Nevertheless, the 
comparison does hint that we should not be surprised to 
find the Ising and Bloch configurations to be competitive 
here. Finally, for the Bloch 180°{110} FDW, which can 
be regarded as one R71° FDW plus one R109° FDW, the 
LDA calculations indicate that the sum of the energies of 
these two walls is much lower than that of a single Ising 
180° FDW. The above caveat has perhaps even more 
force here, but again we can roughly understand in these 
terms why the R180°{110} FDW can only adopt a Bloch 
form. 
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VI. SUMMARY 

In conclusion, we have calculated the domain wall 
widths, energies, and polarization profiles for Ising 
and Bloch ferroelectric domain walls in the zero- 
temperature rhombohedral phase of BaTi03 using both 
first-principles and Ginzburg-Landau-Devonshire meth- 
ods. The first-principles results confirm the expectation 
that 180° domain walls are of Bloch type, adopting a con- 
figuration resembling a pair of 109° and 71° walls in close 
proximity. For the case of the 109° wall, Ising and Bloch 
configurations are competitive. The Ginzburg-Landau- 
Devonshire results are brought into improved agreement 
with the first-principles calculations if the coefficient of 
the gradient term is reduced by about 40%. In view of 
the uncertainties in the original extraction procedure for 
the coefficients, it is not surprising that these parameters 
can be improved; indeed, it is encouraging that even the 
original parameters gave qualitatively sound results. 



While we have not extended our work to other rhombo- 
hedral ferroelectrics such as KNb03, we expect that these 
may show a similar pattern of behavior. When instabil- 
ities other than ferroelectric ones are also present, the 
domain-wall behavior can become more complicated, as 
for example with the octahedral rotations and magnetic 
ordering that play a role in BiFe03. 8 However, we hope 
that the present work will serve as a useful benchmark 
for domain walls in rhombohedral ferroelectrics generally, 
and will lead to an improved understanding of ferroelec- 
tric domain dynamics and switching in these systems. 
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